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We study the infrared behavior of the effective Coulomb potential in lattice SU(3) Yang-Mills 
theory in the Coulomb gauge. We use lattices up to a size of 48 4 and three values of the inverse 
coupling, = 5.8, 6.0 and 6.2. While finite-volume effects are hardly visible in the effective Coulomb 
potential, scaling violations and a strong dependence on the choice of Gribov copy are observed. 
We obtain bounds for the Coulomb string tension that are in agreement with Zwanziger's inequality 
relating the Coulomb string tension to the Wilson string tension. 
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I. INTRODUCTION 

Confinement is one of the peculiar features of Quan- 
tum Chromodynamics, the theory of strong interactions. 
Thanks to 25 years of intensive research in the field of lat- 
tice gauge theory, a few mechanisms for confinement have 
been identified. These mechanisms are associated either 
with monopoles or vortices and seem to be closely related 
to each other In this context, confinement was and is 
stated mostly by a non-vanishing string tension <7wnson 
which expresses the minimal energy of the gluon field 
between a pair of static quarks. The string tension is de- 
fined by Wilson loops and can be extracted in the limit 
of large Euclidean time from the Wilson loop's area-law 
decay. This definition, however, is not completely satisfy- 
ing, because not only quarks but also gluons are confined, 
and there is no area law in the more realistic case of light 
dynamical quarks present in the vacuum. 

There are two other, though less popular, approaches 
that might help to shed additional light on the phe- 
nomenon of confinement. One is given by the Hamil- 
tonian approach which promises to present an under- 
standing not only of bound states but also of the vac- 
uum structure in terms of wave functionals. The other is 
a more field-theoretically inspired approach that focuses 
on the QCD Green's functions and their infrared behav- 
ior. The QCD Green's functions may serve as input to a 
hadron phenomenology based on the Bethe-Salpeter and 
Faddeev equations. There, the ultimate goal is to attain 
a coherent description of hadronic states and processes 
based on the dynamics of confined gluons, ghosts and 
quarks (see, e.g., Ref. 0). 

Both the Hamiltonian approach and investigations of 
QCD Green's functions require to fix a gauge. This in- 



troduces the well-known Gribov ambiguity present in the 
Coulomb as well as in covariant gauges. One should 
keep in mind that the confinement mechanisms associ- 
ated with monopoles and vortices, that received credit 
by reproducing the Wilson string tension, also mostly re- 
quire a gauge condition. In the Coulomb gauge, the Gri- 
bov ambiguity represents a severe source of uncertainty 
and its effect on the results must be faithfully checked. 
On the other hand, the Coulomb gauge yields a partic- 
ularly interesting confinement picture called the Gribov- 
Zwanziger scenario [3J, |J|. This scenario might provide 
an understanding of confinement even in the presence of 
dynamical quarks when the Wilson-loop criterion fails. 

A central element of the Gribov-Zwanziger confine- 
ment scenario in Coulomb gauge is the instantaneous 
color-Coulomb potential involving the Faddeev-Popov 
operator M (in Coulombgauge) and the infrared spectral 
properties of the latter [jj, Q . The expression 

V Coui (x - y)S ab = ( 5 2 [M- 1 (-A)M- 1 ] a6 (x,y)) (1) 

is defined through the vacuum expectation value of the 
potential part of the Hamilton operator 
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resulting from the elimination of longitudinal degrees of 
freedom. Here the potential term Hq u\ is expressed 
in terms of the color charge density (including external 
sources and the charge density of the gluon field itself) 
by means of the color-Coulomb potential, 
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d 6 xd 6 yp a (S)V C oui(x - y)S ab P b (y) . (3) 
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As Zwanziger has shown [|| , the Coulomb potential does 
not equal the Wilson potential Vwiison used to extract the 
string tension awiison as an order parameter for confine- 
ment. Instead, for large spatial distances r the Coulomb 
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potential represents an upper bound for the rise of the 
Wilson potential, 



Vwilson(r) < --VboulM 



(4) 



In other words, there is no confinement without Coulomb 
confinement since the Coulomb stringtension is an upper 
bound for the Wilson string tension [7| , 



^Wilson < gCTCoul 



(5) 



Zwanziger has continuously developed the confinement 
scenario originally proposed by Gribov He has put 
forward the Coulomb potential as a new order parameter 
for confinement 0,0,11]. In fact, the Coulomb potential 
is expected to linearly rise at large r even in the pres- 
ence of dynamical quarks when the Wilson-loop criterion 
fails. Recent lattice studies have shown, however, that 
the relation (|4|) is only a necessary 0] condition for con- 
finement, and that the Coulomb potential can be linearly 
rising with spatial distance even in the dcconfincment 
phase @,[T|. 

Usinglattice techniques, a linearly rising Coulomb po- 
tential [9|, E3, EH and a connection between the center- 
vortex mechanism and the Gribov- Zwanziger scenario 0, 
EH, EH have been observed. Furthermore, Greensite et al. 
proposed [l4[ to use correlators of partial Polyakov loops 
to measure the Coulomb potential. Corresponding SU (2) 
as well as SU (3) studies revealed that the Coulomb string 
tension ercoui could well be 2-3 times lar ger t han the Wil- 
son string tension CTwiison [ToL fill EH, TUl- This is in 
contrast to results of SU{2) studies where the Coulomb 
potential was measured by means of its very definition via 
Eq. (|T} suggesting crcoui = ^Wilson [H, El| • In the present 
study we provide a (yet missing) thorough measurement 
of the Coulomb potential in SU(3) gauge theory based 
on its very definition in Eq. (fT]). We investigate the rela- 
tion between crwiison and ccoui and find, though hedged 
with large uncertainty, crcoui to be 1.6 times larger than 
cwiison- The origin of the systematic uncertainty will be 
discussed. 

The paper is organized as follows. In Sect. [Til we de- 
scribe the details of our numerical simulation and define 
the lattice observables measured. We investigate finite- 
volume effects, lattice-spacing effects and the effects due 
to the Gribov ambiguity in Sect. Hill In Sect. IIVI we 
analyze the infrared behavior of the effective Coulomb 
potential. A summary concludes this paper. 



II. DETAILS OF THE NUMERICAL 
SIMULATION 

A. Lattice samples and gauge-fixing algorithms 



always start from non-gauge-fixed SU (3) gauge configu- 
rations and apply the Coulomb gauge condition subse- 
quently. Our sets of gauge configurations were generated 
with Wilson's one-plaquette action at three values of the 
inverse coupling, f3 — 5.8, 6.0 and 6.2, for a couple of lat- 
tice sizes x L t where L t and L s denote the spatial and 
temporal lattice extension, respectively. We have only 
considered hyper-cubic lattices with L s = L t = L = 12, 
16, 24, 32 and 48. Those ensembles were then gauge- 
fixed to the Coulomb gauge by minimizing the gauge 
functional 



Fu\g] 



(6) 



that involves all space-like links on the lattice. This 
was accomplished by adjusting the gauge transformations 
g x £ SU(3) while keeping the original gauge configura- 
tion U fixed. Due to the particular form of Fu[g] no 
condition is imposed on time-like links. Consequently, 
the different time-slices can be minimized independently. 
We considered gauge-fixing within a given time-slice suc- 
cessful as soon as the stopping criterion 



max Tr 

x, t fix 
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< 10 



-13 



(7) 



was satisfied. Here the lattice gauge-potential is defined 
in the usual way as 



x+i/2,i 
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where a U x ^i = QxUx^g^^ , a is the lattice spacing and go 
the bare coupling constant which is related to (3 through 

To minimize the gauge functional we used an over- 
relaxation (OR) algorithm preceded by an optimally- 
tuned simulated annealing (SA) algorithm. In what fol- 
lows, we call this particular combination of simulated an- 
nealing and over-relaxation steps the SA-OR algorithm. 
To assess the influence of Gribov copies, we also gener- 
ated some gauge copies with the pure OR algorithm with- 
out preconditioning. In all cases, the over-relaxation pa- 
rameter was tuned tow = 1.70 on the small and uj = 1.60 
on the large lattices. More details are given below. 

The SA algorithm has been proven to be very useful in 
handling various optimization problems. For this algo- 
rithm the gauge functional Fi/[g] is regarded as a "spin 
Hamiltonian" where the gauge transformation fields g x 
take the role of "spin variables" coupled through the 
links U x ,i (kept fixed). Minimizing Fu[g] is achieved by 
adiabatically lowering the auxiliary temperature T of a 
statistical spin glass system characterized by the Gibbs 
weight 



W\g]<xexp(-Fu\g]m 



(9) 



For our study we use the standard lattice formulation 
of SU(3) Yang-Mills theory in Coulomb gauge where we 



The minimization process always starts with equilibrat- 
ing this spin system at some initial temperature T = Ti 
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which is then slowly decreased. Formally, in the limit of 
(adiabatically) lowering T — > this system approaches 
the ground state and hence the gauge functional reaches 
its absolute minimum for a given gauge configuration. 
For the practical purpose considered here such an adi- 
abatic cooling-down process is not feasible as it would 
require an enormous amount of computing time. Never- 
theless, we find that much lower minima for i*i/[</] can 
be reached, compared to applying only over-relaxation 
(OR), if we combine the SA with the OR algorithm as 
follows: We start from an initial temperature of T = 0.45 
and linearly decrease the temperature down to Tf = 0.01 
within 1500 "compound sweeps" . Each such sweep con- 
sists of one heatbath and three microcanonical update 
sweeps. After this, we use the OR algorithm until the 
Coulomb gauge is reached, i.e. the stopping criterion ([7]) 
is satisfied. 

The advantage of using the SA-OR instead of the OR 
algorithm alone becomes more pronounced as the lattice 
becomes larger. Furthermore, the number of necessary 
iterations in the subsequent OR algorithm is drastically 
reduced by a preceding SA algorithm, the more the lower 
the final Tf is chosen. Note that instead of adding sub- 
sequent OR steps, we could also have used SA on its own 
extending its use to a much lower temperature Tf to fix to 
Coulomb gauge. This, however, is much more CPU-time 
intensive and we find no benefit in doing this, because 
after gauge-fixing the transversality condition ([7]) must 
at any case be guaranteed with high precision, which can 
be achieved only by the finalizing OR. 

We observe that the time-slices of a given configuration 
may behave very differently during the iterative gauge- 
fixing process. In fact, we find the number of necessary 
iterations may differ by a factor of 10 to 20 between the 
individual time-slices of a given configuration. In the ma- 
jority of cases, time-slices did not show any recalcitrancy 
during gauge fixing, although in some cases time-slices 
could not be fixed within a certain (predefined) number 
of iterations. In the latter case we simply repeated the 
entire gauge-fixing process for these time-slices, using the 
same algorithm but starting from a different randomly 
chosen gauge transformation. The "well-behaved" and 
hence already gauge-fixed time-slices were not touched 
again. 

After all individual time-slices had been minimized, 
the original configuration U was gauge-transformed, i.e., 



g u x 



To simplify the notation we drop the 



label g in what follows and assume that a gauge configu- 
ration U satisfies the Coulomb gauge condition already. 
Since our observables, namely the effective Coulomb po- 
tential and the ghost propagator, are genuine three- 
dimensional, instantaneous observables defined by space- 
like links only, we did not have to fix the residual gauge 
freedom. The latter, after the Coulomb gauge has been 
fixed, resides in spatially constant but time-dependent 
gauge transformations (for a continuum view at this 
problem see [17|). 



B. Observables of interest 

The Coulomb energy is a complicated functional of 
the transverse gauge potential Ai(x) and the total color 
charge density. Nevertheless, it is instructive to charac- 
terize its gross features through the infrared and ultravio- 
let behavior of the expectation value of the color-diagonal 
part of the kernel M _1 (— A)M _1 in momentum space 
alone. On the lattice this is defined as the MC average 



Coul 



(k) 




ik-(x-y) 



M- 1 {-A)M- 



xy 



(10) 

where we use a shorthand notation for the scalar prod- 
uct k ■ x — 2ir X}j=i kiXi/Li with integer- valued lattice 
momenta fcj and lattice coordinates Xi- M is the lattice 
Faddeev-Popov operator for the Coulomb gauge 

3 

M£ = S Xi , Vi ]T <He Tr [ {T a , T b } (u x>i + U x _^ 6g,g 

2=1 



2 TT" U x i - - 2 T a T U 



(ii) 



Note that the Faddeev-Popov operator is a direct sum 
of operators acting within individual time-slices. In co- 
ordinate space these three-dimensional operators define 
the Coulomb energy of a given dynamical (gluonic) color- 
charge density plus an external one (cf. Eq. ([3])). Given 
the tree-level form of the Coulomb potential on a three- 
dimensional lattice we relate integer-valued lattice mo- 
menta ki S (—Li/2, Li/2] to physical ones by 



qi{ki) = - sm — 
a V Li 



(12) 



Physical units are assigned by using the interpolation 
formula for r^/a as given in [l8| setting r§ = 0.5 fin. To 
simplify the writing we introduce q as abbreviation for 
\q\ whenever appropriate. 

In Ref. [H an analytic calculation of the Coulomb po- 
tential is presented which reads, upon Fourier transfor- 
mation, 



V G oul(q)=q 2 G 2 (q)+V c (q) 



(13) 



Here G denotes the ghost propagator (entering the dis- 
connected part) and V c denotes the connected part of 
the potential. Under the assumption that the (yet un- 
known) connected part can be neglected, an infrared 
asymptotic limit for Vboul has been given in [gj. It will 
be analyzed below at what momenta the factorization 
Vboui(<7) — q 2 G 2 {q) is justified from our data concerning 
both the effective Coulomb potential and the ghost prop- 
agator. The latter can be estimated in momentum space 
as the MC average 



G L (k) 




a ik-(x-y) 



[M 



aa 
J xy 



(14) 
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at non-zero lattice momenta k. As for the Coulomb po- 
tential we use Eq. (TT2]) to assign physical momenta to 
G. To invert the Faddeev-Popov operator we adapted 
the techniques developed in Landau-gauge studies of the 
ghost propagator (see, e.g., [Hj]). The data for the ghost 
propagator used to test the factorization hypothesis will 
not be presented in the present publication. They have 
been presented at Lattice 2007 [20J and will be discussed 
more in depth in a forthcoming paper (2lj . 

Note that both the evaluation of the effective Coulomb 
potential and of the ghost propagator involve CPU-time 
intensive operations. As a consequence, we have re- 
stricted our measurements to lattice momenta k that sur- 
vive a cylinder cut. Our cylinder cut is the obvious adap- 
tation of the Landau-gauge cylinder cut [22^ . To mini- 
mize finite- volumes effects, we also cone cut our data [22] 
if they refer to lattices smaller than (2.5 fm) 4 . 

C. Running coupling and physical scale 



accessible in our simulations really permit a feasible fit to 
the one- loop or the two-loop expression given in Eq. (fTSj) . 
We find that this is not the case and that the ultraviolet 
fit described in [2(| has artificially up-scaled our data by 
a free factor bigger than one. In the present study, we 
therefore do not rely anymore on this ultraviolet fit. 

Indeed, the physical scale is fixed by simply multiply- 
ing the bare lattice data for the effective Coulomb poten- 
tial with 6/ (/3a 2 ) 

^Coui(o) = |a 2 V c L oul (fc,/3) , (19) 

where a denotes the lattice spacing in GeV" 1 . A gam, we 
use the interpolation formula in [18| to set a assuming 
r = 0.5 fm. For all figures in the present paper, the 
physical scale of the effective Coulomb potential is fixed 
in this way. 



The Coulomb potential is a renormalization-group in- 
variant which can be written as (here for pure SU(3) 
gauge theory) [23| 



o 2 Vcoui(o) 



12 
IT 



ffCoulteMcoul) 



(15) 



where Acoui is a special QCD scale parameter charac- 
teristic of the Coulomb gauge, that defines a running 
coupling constant <?Coui- The latter has to satisfy the 
renormalization-group equation 



1 — -Q- — — PCouliffCoulJ , 



(16) 



where the beta function, /?coui, has the usual weak- 
coupling expansion starting with the two standard 
scheme-independent coefficients 



b 



11 
16tt 2 



and bi = 



51 



1287T 4 



(17) 



(see [23| for higher terms). For sufficiently large q, 
the product llg 2 Vc ui(g)/12 is expected to be described 
through the two-loop expression of the running coupling 



2 

ffCoul 



(<?) = 



6 ln(o 2 /A 2 oul ) 



0l ln[ln(g 2 /A 2 oul )] 



^5 



ln(q 2 /A 2 



CouV 



Lattice data describing o 2 Vcoui(o) do not depend on 
the lattice spacing a in the asymptotic scaling region. At 
larger a, scaling violations should be expected though, 
and they will be discussed below for the lattice spacings 
used in this study. 

In a previous analysis of the data [2(| we used an ul- 
traviolet fit to the one-loop expression (cf. the first term 
of Eq. (fT8")0 to fix the unknown physical scale of the ef- 
fective Coulomb potential (see Ref. [2(| for details). For 
the present study, we scrutinized if the highest momenta 



III. STUDYING SYSTEMATIC EFFECTS 

In this section we discuss the effects of finite lattice 
volumes and lattice spacings as well as the influence of 
the Gribov ambiguity on the effective Coulomb potential. 



A. Lattice artifacts 

As we are primarily interested in the product 
o 4 Ucoui(o), we directly discuss this product instead of 
the effective Coulomb potential itself. Note that we in- 
vestigate effects of finite lattice spacings and volumes 
by considering Coulomb-potential data collected for first 
SA-OR copies only. 

Finite- volume effects are studied by varying the lattice 
sizes from 12 4 to 48 4 but keeping /3, and hence the lat- 
tice spacing a, fixed. We find that only data obtained on 
the smaller lattices, 12 4 and 16 4 , at (3 = 6.0 show visible 
finite-volume effects at lower momenta. For larger lat- 
tices, namely 24 4 , 32 4 and 48 4 , effects seem to be mild 
(see the magnified view in the right panel of Fig. [1]). At 
(3 = 5.8 (reaching even lower momenta) only the 12 4 data 
clearly deviate from the other data (see Fig.[TJ left panel). 

Lattice-spacing effects are investigated by comparing 
data from lattices of equal physical volume for different 
values of [3. Within our choice of (3 values and lattice 
sizes L 4 , we can find only a few combinations of f3 and L 
where this is approximately possible. For those we can 
compare data at approximately equal physical momenta 
and disentangle by eye the effect of varying a. As demon- 
strated in Fig.[5]these discretization effects are small and 
of the order of 10 to 15% at largest. The difference be- 
tween data for (3 = 6.2 and 6.0 is smaller (right panel) 
than the difference between data for (3 = 6.0 and 5.8 (left 
panel) . 
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FIG. 1: (Color online) The Coulomb potential multiplied by q 4 shown as a function of q 2 in physical units. We show data for 
different lattice sizes at /3 = 5.8 (left) and f3 = 6.0 (right) to illustrate finite-volume effects. These seem to be under control for 
data on lattices larger than 16 4 because those fall roughly on the same curve. For both j3 values we only used data from first 
SA-OR copies (fc). 
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FIG. 2: (Color online) The Coulomb potential multiplied by q 4 shown versus q 2 measured on comparable physical volumes. 
Data for /3 — 5.8 and 6.0 are shown on the left, and for /3 = 6.0 and 6.2 on the right hand side. Only data from first SA-OR 
copies is shown. 



B. Effects due to the Gribov ambiguity 

In comparison to lattice artifacts, the Gribov ambi- 
guity turns out to have a much larger impact on the 
Coulomb potential data. In order to assess the influence 
of this ambiguity we follow the "first copy - best copy" 
(fc-bc) strategy applied before in Landau-gauge studies 
of gluon and ghost propagators [H, H3|- Here, we use 
this strategy in two different ways. 



First, we estimate the number N cp of gauge-fixed 
copies per configuration necessary to achieve "quasi con- 
vergence" of the Coulomb potential, considered as a func- 
tion of iV C p. Second, we quantify the systematic error of 
the Coulomb potential that is admitted if an arbitrary 
(first) gauge-fixed copy is chosen instead of the best copy 
among N cp copies for each gauge configuration. A copy 
is considered to be the best among all N cp gauge-fixed 
copies of a given configuration if its gauge functional is 
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FIG. 3: (Color online) The upper row shows the convergence of Vcoui with increasing number N cp of copies obtained with 
either the SA-OR (circles) or the OR algorithm (crosses). The Coulomb potential is shown for the smallest (on-axis) lattice 
momentum available for each of the three lattice sizes, i.e. k — ([1,0,0]) with square brackets indicating that the average 
over all permutations is taken. The lower row shows the corresponding convergence of the (average) deviation of the gauge 
functional from its best value known for each configuration. As a function of A'cp < A™ ax , the values of Vc ou i and F are to 
be understood as those for the "currently best" gauge copy among the iV cp copies inspected for each configuration after N cp 
repetitions of gauge fixing. 
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FIG. 4: (Color online) The ratio of the Coulomb potential (as function of q 2 ) if evaluated either on arbitrary (first) copies from 
simple OR or on best gauge copies from repeated SA-OR. Note that the enhancement for the first copy tremendously grows 
with q 2 — *■ compared to what is known in the case of the ghost propagator. In contrast to the upper panel of Fig. [3] only 
momenta are included that are allowed by the cuts mentioned above. 
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lower than those of all the other N cp — 1 copies after the 
gauge-fixing has been attempted iV cp -times. 

Let us first compare the convergence of the bare data 
describing the Coulomb potential upon increasing N cp — > 
7V™ ax for the two gauge-fixing algorithms OR and SA- 
OR. As an example, in Fig. [3] we show for each lattice 
size data describing Vboul f° r the lowest (on-axis) lat- 
tice momentum available, i.e. k = ([1,0,0]) with square 
brackets indicating that the average over all permuta- 
tions is taken. Note that in contrast to the rest of this 
paper we did not apply neither the cylinder nor the cone 
cut here. The obtained deviations from the best-copy 
value can be considered as an upper bound for all other 
momenta. 

The data were obtained at (3 = 6.0 on lattice sizes 
L 4 = 12 4 , 16 4 and 24 4 . For each gauge-field configura- 
tion a number of A™ ax = 20, 30 and 40 independent 
gauge-fixed copies was generated separately with the OR 
and the SA-OR algorithm. From the figure we see that 
upon increasing N cp the effective Coulomb potential de- 
creases and becomes (more or less) independent of iV C p 
for A cp coming closer to iV™ ax . What A cp is sufficient to 
achieve "quasi convergence" depends, of course, on the 
gauge-fixing algorithm and on the lattice size. In fact, 
it is clearly visible in Fig. [3] that the number of gauge 
copies necessary to achieve convergence is substantially 
lower for the SA-OR algorithm than for the OR algo- 
rithm. With both algorithms one needs to consider more 
gauge copies with growing L. For example, if we use the 
SA-OR algorithm, a number of copies A cp = 5 on a 12 4 
lattice and N cp = 15 on a 24 4 lattice is sufficient. In 
contrast, if we were using the OR algorithm, A cp = 40 
or more copies are necessary for a lattice like 24 4 . Note 
that the observed increase with L at fixed (3 is partly due 
to the smaller physical value associated with the lowest 
on-axis lattice momentum that needs to be considered 
with increasing L. 

On the larger lattices 32 4 and 48 4 we could not afford 
to gauge-fix more than a single gauge copy per configura- 
tion with the SA-OR algorithm. This was simply due to 
a drastic increase of the necessary number of iterations, 
but also due to a growing number of "trouble-making" 
time-slices encountered on those larger lattices. Thus, 
we did not apply the OR algorithm for the purpose of 
comparison, and therefore we are not in the position to 
assess the influence of Gribov copies on these lattices at 
the present stage. 

As mentioned above, we also used the fc-bc strategy 
to estimate the impact of Gribov-copy effects on the 
Coulomb potential data at different physical momenta. 
For this purpose we gauge-fixed our field configurations 
at (3 = 5.8, 6.0 and 6.2 to Coulomb gauge only once with 
the OR algorithm on one hand and 7V cp -times with the 
SA-OR algorithm on the other. In order to obtain the 
results shown in Fig. [4] we have chosen, refering to Fig. [3l 
A cp = 10, 15 and 20 as sufficient numbers of gauge-fixed 
copies per configuration for the lattice sizes 12 4 , 16 4 and 
24 4 , respectively. Then the Coulomb potential was mea- 



sured separately on the set of single copies obtained with 
the OR algorithm on one hand and on the set of best 
copies obtained with the repeated SA-OR algorithm on 
the other. For brevity we refer below to these two sets 
as the first OR and the best SA-OR copies. 

The ratio of the effective Coulomb potential measured 
for first OR copies and best SA-OR copies is depicted 
in Fig. |H as a function of momentum squared. The 
Gribov ambiguity has a dramatic impact on the effec- 
tive Coulomb potential at q 2 < 10 GeV 2 . Even for the 
rather small lattices considered here, and hence for rather 
high physical momenta, the measurement of the effective 
Coulomb potential on first OR-copies gives results larger 
by up to 100% than the results on best SA-OR copies. 
Note that this effect is much stronger than what has been 
observed for the ghost propagator using the same method 
[20l |2lj . There an enhancement of about 5 to 10% was 
typical for the presently accessible lowest momenta. Note 
also that here we are comparing the standard with one 
of the best presently known methods of gauge fixing. 

In order to assess next the difference between Vcoui ob- 
tained from an arbitrary, first SA-OR copy (fc) and the 
best among a sufficient number of SA-OR copies (be), 
restricted, however, to the smallest lattice momentum 
for each lattice size, i.e. k — ([1,0,0]), we have to look 
back to Fig. [3J Considering the ratio R of the effective 
Coulomb potential measured either for first or best SA- 
OR copies as a function of the lattice size L, we find 
that this is well described by R ~ c — d/L. If such an 
ansatz was used to extrapolate the ratio R at (3 = 6.0 
to L = 48, a ratio R = 1.6 would be obtained. For the 
first OR copy an overestimation factor R = 2.7 would be 
expected. Both are reasonable upper bounds for the over- 
estimation of the Coulomb potential at any fixed physi- 
cal momentum for first SA-OR copies and - even worse - 
first OR copies. This estimate will be needed in the next 
Sect. [TV] 

We conclude that the effective Coulomb potential is 
less affected by Gribov-copy effects if we use the SA-OR 
algorithm instead of the OR algorithm. This conclusion 
rests on the observation that the "quasi convergence" of 
the Coulomb potential is faster for the SA-OR than for 
the OR algorithm. Second, the results obtained on ar- 
bitrary, first SA-OR copies (fc) are less affected by the 
Gribov ambiguity than those obtained on arbitrary, first 
OR-copies. Therefore, we have used the SA-OR algo- 
rithm as our method of choice for the results to be pre- 
sented in the following. Recall, however, that if only first 
SA-OR copies are available for analysis, measurements of 
the Coulomb potential in the infrared region will be ac- 
companied with an increased uncertainty. For instance, 
an overestimation of about 60% for the smallest lattice 
momentum on a 48 lattice at (3 = 6.0 must be expected. 
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IV. INFRARED BEHAVIOR 

Despite the Gribov ambiguity being that large, we now 
try to summarize what we know about the momentum 
dependence of the effective Coulomb potential, globally 
and in particular in the low-momentum region. As men- 
tioned above, we were not in the position to generate 
more than a single SA-OR gauge copy per configuration 
on the larger lattices 32 4 and 48 4 . Therefore we present 
here a full set of data concerning the Coulomb potential 
for a single SA-OR copy (fc) per configuration for all (3 
values and lattices sizes, ensuring in this way an equal 
treatment of Gribov copy effects on both small and large 
lattices. As is well known, this choice is equivalent to 
an averaging over all local minima of all configurations, 
i.e. all over the Gribov Region. Best-copy data, that 
we have available only up to lattices 24 4 (after inspect- 
ing a sufficient number N cp of copies) would come close 
to a prescription that requires an average over only the 
absolute minimum per configuration, i.e. restricted to 
the Fundamental Modular Region. Zwanziger has argued 
that these averages should approach each other in the 
limit of large volumes. 

As long as they did not converge, we are admitting a 
strong systematic effect when we restrict the analysis to 
the first SA-OR copy. This can be clearly seen in our 
data from smaller lattices, 12 4 , 16 4 and 24 4 . In order not 
to overload the Fig. [5] we show here (and lateron in Fig. [6]) 
only selected results from SA-OR best-copies. The data 
are from the 24 4 lattice where we had the choice between 
40 copies at (3 = 6.0 and between 20 copies at [3 = 5.8. 
We try here (and later for the Coulomb potential) our 
best to estimate the systematic error emerging from the 
ignorance of further Gribov copies at larger lattices. A 
thorough study of the Gribov ambiguity for larger lattices 
remains highly desirable. 

In the infrared momentum region, the running cou- 
pling given through the effective Coulomb potential di- 
verges stronger than l/q 2 . This is shown in the left panel 
of Fig. [5l The very fact of the infrared enhancement will 
not need to be revised if once a systematic account for 
the Gribov effect would be undertaken, although the di- 
vergence would be less pronounced. A rough indication 
of the size of the effect is given by the filled symbols in 
that figure. These are best-copy results for the lattice 
24 4 . 

In Ref . [8] , Zwanziger presented an analytic calculation 
of the Coulomb potential. By only considering the dis- 
connected part of the expectation value of the effective 
Coulomb potential (cf. Eq. Q2])) Zwanziger predicted an 
almost linearly rising effective Coulomb potential in the 
infrared limit. Using our data for the ghost propaga- 
tor [2(| we are now in the position to test the validity of 
his factorization hypothesis. If Zwanziger's assumption 
were valid, the ratio 



should be constant as function of the momenta. Note 
that V L and G L denote the bare lattice Coulomb poten- 
tial and the lattice ghost propagator taken at the physi- 
cal momentum q. The resulting plot shown in the right 
panel of Fig. [5] demonstrates that the assumption of fac- 
torization is valid only for q 2 > 10 GeV 2 , but it is not 
correct in the momentum range q 2 < 10 GeV 2 . This is 
in agreement with the results of Langfeld and Moyaerts 
for SU(2) pure gauge theory [l6[. 

Much alike the enhancement of the running coupling, 
our conclusion that the factorization hypothesis is vio- 
lated would also not be invalidated if the effect of Gribov 
copies was taken into account properly. Similar to the 
left panel, the anticipated Gribov effect on the violation 
of factorization is shown by the filled symbols in the right 
panel, which are the best-copy results for the lattice 24 4 . 

We discuss now the momentum dependence of the ef- 
fective Coulomb potential in three clearly emerging mo- 
mentum ranges, the high-momentum range, the interme- 
diate momentum range and the low-momentum range, 
and describe the influence of the Gribov ambiguity in 
each range separatly. 

The left panel of Fig. [S] shows the presently 
known picture concerning the momentum dependence of 
q 4 Vcoui(q)- A logarithmic momentum scale has been 
chosen in order to give a global view including the ultra- 
violet and infrared behavior. For the largest momenta in 
the high-momentum range q 2 > 10 GeV 2 , the Coulomb 
potential shows roughly the expected l/q 2 behavior lead- 
ing to an increase of q 4 Vc ui(<7) linear in q 2 . From Fig.0] 
it is clear that the high-momentum region is robust with 
respect to the Gribov ambiguity. Although the inspec- 
tion by eye suggests that we are seeing the tree-level form 
of the Coulomb potential, we could not find reasonable 
fits of our data by the one-loop or the two-loop expres- 
sions given in Eq. (fT8)) . We conclude that much higher 
momenta must be considered to get an estimate of the 
Coulomb scale parameter Ac ou i from such a fit. 

With decreasing physical momenta, the first-copy data 
for q A Vcoui('z) reach an almost flat region in the inter- 
mediate momentum range 0.2 GeV 2 < q 2 < 6 GeV 2 , al- 
though a little bulge is visible in the left panel of Fig. [5] 
If the function q 4 Vcoui(<?) stayed constant on the level 
of 20 GeV 2 in the limit q 2 — > 0, this would imply a 
perfect linearly confining potential corresponding to an 
estimate of acoui ~ (890 MeV) 2 . This figure is more 
likely an upper bound. 

Indeed, if for large spatial distances r we assume the 
simple ansatz [15| 



Vboui(r) = ~cr Coul r + C/r , 
this suggests a momentum behavior 

87T(7Coul , 47rC 



VcouiO?) 



(21) 



(22) 



(aqG L (q)) 2 



(20) 



with the intercept of q 4 Vc u\(q) at g 2 = defining the 
Coulomb string tension ccoui- 
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FIG. 5: (Color online) Left: The running coupling f/c ou i(g) K 9 2 ^Coui(q) diverges in the infrared region and tends to zero in the 
asymptotic limit q 2 — * oo. Right: The factorization of the effective Coulomb potential is violated for momenta 0.04 GeV 2 < 
q 2 < 10 GeV 2 . For both the running coupling and the test of the factorization hypothesis, the open symbols (including stars) 
represent measurements on the first SA-OR copies per configuration. For comparison, the filled triangles and filled squares 
show selected results for the best SA-OR copy (be) per configuration for two (3 values on the largest lattice 24 4 where the 
Gribov problem was fully under control. 



In the intermediate momentum range 0.2 GeV 2 < q 2 < 
6 GeV 2 the Gribov effect sets in and becomes appar- 
ently more severe with decreasing momentum. For in- 
stance, for the smallest (on-axis) lattice momenta on lat- 
tices of sizes 12 4 , 16 4 and 24 4 we have seen in Fig. [3] 
that the Coulomb potential Vc ui(<7) is overestimated by 
the first SA-OR copies compared with the best SA-OR 
copies (among 40 copies). For the 24 4 lattice at (3 = 6.0 
the overestimation amounts to » 40%. This can be ex- 
trapolated to the 48 4 lattice where the effect amounts 
to « 60%. This is an upper bound for the Gribov ef- 
fect experienced by Vc ul(?) at physical momenta that 
are allowed by the cylinder and cone cuts. 

In agreement with these estimates it can be seen in the 
left panel of Fig. [5] that in the intermediate momentum 
range 0.2 GeV 2 < q 2 < 6 GeV 2 the best-copy data from 
the 24 4 lattice (shown as filled symbols) provide us with 
another, independent early indication of a plateau. The 
somewhat lower level of « 10 GeV 2 would correspond to 
CTCoui ~ (630 MeV) 2 . In view of this the bulge must be 
understood as an artefact of insufficient gauge fixing. 

With the simulation reported here, on our largest lat- 
tices the low-momentum range with q 2 < 0.2 has become 
accessible for the first time. Rather unexpectedly in this 
region the first-copy data for q 4 Vcoui(q) drop with de- 
creasing momentum as seen in the left panel. The right 
panel of Fig.[B]shows the infrared region magnified and in 
a linear scale in q 2 . This picture shows that a fit ansatz 
linear in q 2 describes the drop of the first-copy data very 
well. We do not know whether a similar effect, namely 
the onset of an apparently new infrared regime in the 



low-momentum range will happen for the best-copy data 
as well. For the time being we assume that the bulge and 
the new infrared regime is only a matter of measurements 
on insufficiently gauge-fixed configurations. We have fit- 
ted the behavior according to Eq. (f2"2")l . The Coulomb 
string tension is estimated as 

cr Cou i = (552 ± 35 MeV) 2 . (23) 

With some caution we may consider this as the common 
limit for q 2 — > and the common lower bound for the 
Coulomb string tension (common to both standards of 
gauge-fixing). 

The other fit parameter, the "Coulombic" coefficient 
C in front of the 1/r term in Eq. (|2ip . is obtained as 

C = 6.0 ±1.0. (24) 

This parameter has no relation to the "Coulombic" part 
1/r in Eq. (|21|) . It rather describes the narrow momen- 
tum interval where the single-copy data probably con- 
verge to the best-copy results for q 4 Vboul(?)- The small 
number of data points is another reason why we give not 
much significance to the fit. Still, details of the least \ 2 
fit of the first-copy data are shown in Table HI We re- 
mark that the choice of the upper momentum cutoff for 
the fitting range, q max , has only weak influence on the fit 
results. 

In units of the Wilson string tension the fit result is 
ccoui = (1-6 ± 0.2) awuson ■ (25) 

This is the tentative lower bound for the Coulomb string 
tension. 
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TABLE I: Results of X 2 fits to the single-copy data at mo- 
menta q 2 < 5^ ax . 



?L [GeV 2 ] 


# data 
points 


V^Coui [MeVJ 


C 


X 2 /ndf 


0.11 


5 


534(16) 


6.6(3) 


2.9 


0.16 


6 


526(18) 


6.8(4) 


1.8 


0.17 


7 


558(20) 


5.8(2) 


2.5 


0.18 


8 


587(28) 


4.9(2) 


3.8 



Our estimate for the Coulomb string tension is in 
agreement with Zwanziger's inequality. The relevance 
of this estimate is, however, faced with three sources of 
uncertainty. 

• First, it relies only on first-copy data in a rather 
small number of data points, and the obtained 
X 2 /ndf values are rather large (see Table |TJ) . The 
latter might be interpreted as a probable inade- 
quacy of the assumed infrared ansatz Eq. (|22p. On 
the other hand, the right panel of Fig. [5] supports 
such a behavior. 

• Second, the weak but visible scaling violation of 
the effective Coulomb potential has the effect that 
our estimate of the Coulomb string tension would 
be higher if we considered higher inverse coupling 
constants (3. The effective Coulomb potential in 
general slightly increases with increasing 0. 

• Third, the strong Gribov effect is neglected in this 
estimate for the Coulomb string tension. If we had 
consequently looked for the best SA-OR copies and 
had measured q Vcoui(q) for these, the amount of 
overestimation by the first-copy data in the bulge 
region is of the estimated order. One possibility is 
that by the drop described by the fit given above 
the (yet unknown) level of the best-copy results is 
reached. However, we cannot exclude the possibil- 
ity that the best-copy data in the low-momentum 
range will also enter a new infrared regime with a 
similar decrease, such that the overestimation by 
the first-copy data remains. In this case the final 
estimate of the Coulomb string tension would be 
close to the Wilson string tension. 

In the light of these uncertainties, we find it difficult 
to draw a conclusion on the exact value of the Coulomb 
string tension. Our value is larger than the values re- 
ported in previous SU(2) investigations starting also 
from the definition Eq. (fTJ). These authors arrived at an 
estimate close to the Wilson string tension [l5|, [n| . How- 
ever, the Gribov copy problem for the effective Coulomb 
potential was ignored in these studies. Furthermore, in 
Ref. [HI the estimate of the Coulomb string tension ac- 
tually relies on data in the perturbative region, while 
the first plateau of q 4 Vcoui(<z) nas been considered as a 
finite-volume effect. Such a plateau could be observed 



in Ref. [l6| but the further decrease of q 4 Vboul for even 
lower momenta was beyond the possibilities of this in- 
vestigation. In contrast, SU(S) studies using incomplete 
(partial- length) Polyakov lines, made in order to interpo- 
late between the Coulomb string tension and the Wilson 
string tension, gave ctcouI = {2~3)a W iison [19, [HI- These 
studies also have neglected the problem of Gribov copies 
that might have affected the measured correlators. 



V. CONCLUSIONS 

In this study we have attempted a thorough measure- 
ment of the effective Coulomb potential in SU(3) lattice 
gauge theory. We used a broad range of lattice sizes, 
12 4 — 48 4 , to perform Monte-Carlo simulations at the 
three values f3 = 5.8, 6.0 and 6.2. This has allowed us 
to show that finite-volume effects are hardly visible on 
the larger lattices and discretization effects are modest. 
Additionally, the use of the fc-bc strategy has revealed 
a dramatic dependence of the Coulomb potential on the 
choice of Gribov copy. 

Unfortunately, by computer resources we were forced 
to restrict this "Gribov analysis" to the smaller lattices 
12 4 , 16 4 and 24 4 . Thus, performing a full Gribov study 
up to the largest lattices still remains a highly desirable 
goal. We note that the necessity of choosing best copies 
versus first (and hence arbitrary) g aug e- fixed copies is a 
matter of current debate (see also [25l. [26|). As another 
example, in a BRST formulation, an average over all Gri- 
bov copies is taken which, on the lattice, usually leads to 
the well-known Neuberger 0/0 problem [27ll2a|. For a re- 
cent lattice BRST formulation without this complication 
see [29j . 

What can be said here with confidence is that for the 
effective Coulomb potential we find an extraordinarily 
strong Gribov-copy effect which has never been observed 
before for other observables (say, the gluon and ghost 
propagators). 

We see a strong violation of the factorization hypoth- 
esis for the effective Coulomb potential in momentum 
space below q 2 ~ 5 GeV 2 . For smaller momenta the 
"connected part" of the corresponding expectation value 
in Eq. (|13p is not negligible anymore. This spoils any 
simple relation between the momentum dependence of 
the effective Coulomb potential and the behavior of the 
ghost propagator. 

Using only one SA-OR gauge copy per configuration 
and hence allowing a strong systematic Gribov effect to 
be included in the bargain, we found a new infrared 
regime of the effective Coulomb potential. The first 
plateau of q Vcoui(<?)) encountered with decreasing mo- 
menta, turns out not to represent the asymptotic behav- 
ior, because there is a further step-like decline for even 
smaller momenta. The size of the step is of the same 
order as the extrapolated difference between single-copy 
and many-copy SA-OR results. Therefore, we adopted 
the point of view that the "breakthrough" to some new 
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FIG. 6: (Color online) The effective Coulomb potential multiplied by q 4 as a function of the physical momentum squared. 
Left: with a logarithmic momentum scale in order to overlook both the IR and UV behavior. Right: the infrared momentum 
region is shown in a linear scale in q 2 in order to judge the adequacy of the linear fit of the first-copy data in the extremely IR 
region. The infrared fit used to extract the corresponding Coulomb string tension is also shown in the left panel. Open symbols 
(including stars) represent measurements on the first SA-OR copies per configuration. For comparison, the filled triangles and 
filled squares in both panels show results for the best SA-OR gauge copies (be) for two [3 values on the lattice 24 4 , the largest 
lattice where the Gribov problem was under scrutiny. 



infrared regime at large enough volume is a feature only 
of the single-copy data, and that some kind of conver- 
gence (between averaging over the Gribov Region and 
the Fundamental Modular Region) is behind this obser- 
vation. 

Future studies shall scrutinize whether the presumed 
common infrared limit of these two schemes really ex- 
ists or, alternatively, the Gribov ambiguity persists at 
lower momentum for larger volumes. We estimated the 
Coulomb string tension by fitting the data at the lowest 
momenta and found it approximately 1.6 times larger 
than the Wilson string tension. If the Gribov ambiguity 
persists, it is not excluded that in the - further delayed 

infrared limit finally ccoui = o"wiison will be found. 
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